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Abstract. The possibility of a neutron m = 2-superfluid 
in the interior of neutron stars is investigated. This pair- 
ing state is energetically favoured in strong magnetic fields 
(H ~ 10 16 — 10 17 G). Because of the node in the angular- 
dependent energy gap along the field direction the neu- 
trino emissivity is only suppressed polynomially as func- 
tion in T/T c instead of exponentially, as it is obtained 
for a nodeless pairing state. The effect of this pairing 
state on the thermal evolution of neutron stars is stud- 
ied, and its outcome is compared with the evolution of 
"normal" , i.e. nodeless, superfluid and non-superfluid neu- 
tron stars, and also with observations. We find that par- 
ticularly the predicted surface temperatures of the en- 
hanced cooling scenario considerably change and come 
into agreement with temperatures deduced from observa- 
tional data within the hydrogen atmosphere model. Fur- 
thermore the surface temperature depends on the mag- 
netic field strength as an additional parameter aside from 



see also Voskresenskii fc Senatorov, 1986 ) the uncertain- 



the neutron star mass. The latter is however only opera- 
tive in the case of the intermediate cooling scenario. 

Key words: Stars: neutron - Stars: evolution - Dense 
matter - Magnetic fields - X-rays: stars 



1. Introduction 

Comparison of soft X-ray and extreme UV observations 
with simulations of the thermal evolution of neutron stars 
gives us a powerful tool for investigating the interior 



ties in the degree of modifications of the modified Urea 
rate by medium effects suggest a third possible scenario, 
the intermediate cooling scenario. 

It is generally believed that neutrons in the interior of 
a neutron star form superfluid pairs in a 3 P2-state be- 
low some critical temperature T c ~ 10 9 K. Since this 
high density (p > p = 2.8 x 10 14 gcm~ 3 ) matter in (3- 
equilibrium is not directly accessible by laboratory ex- 
periments, one can only extrapolate the two-body inter- 
action to this regime. The pairing gaps however depend 
very sensitively on the underlying microscopic model for 
the effective nucleon-nucleon interaction. Its value and 
the density range over which pairing occurs is therefore 
quite uncertain. Another uncertainty concerning the 3 P2- 
pairing of neutrons refers to the angular dependency of the 
gap. In isotropic matter the m = 0-angular dependency 
A 2 (fc) oc 1/3 + cos 2 9 is energetically favoured (Hoffberg 
Above some critical magnetic field strength 
10 17 G however the free energy is minimised 



et al. 



, 1970) 



by a m = 2-dependency A 2 (fc) oc sin 2 9, where 9 is the an- 
gle between the quasiparticlc momentum k and the mag- 
netic field direction (Muzikar et al., 1980). This differ- 



ent angular dependency turns out to be important for the 
cooling behaviour of superfluid neutron stars, since the 
m = 2-angular dependency has nodes at 8 = 0, tt. The 
available phase space of the various processes keeps there- 
fore undiminished in the direction of the magnetic axis, 
and the reaction rates a re only suppressed polynomially 



such s tars, since young neutron stars cool mainly via nou v v ' 3' 



j instead of exponentially ( Anderson fc Morel, 1961 ; Levcn- 



trino emission from there. The rate of neutrino emission 
depends on a number of ingredients, among which the 
composition and the superfluid pairing play the most im- 
portant roles. 

Differences in the composition of various models make 
it sensible to distinguish between slow (or standard) and 



The observed surface magnetic fields of neutron stars 
ranges from B = 1.7 x 10 s G (PSR B1957+20) up to 
2.1 x 10 13 G (PSR B0154+61), with a median value of 
about 10 12 G. These values are deduced from the observed 
spin-down velocity of 558 pulsars by assuming the mag- 



nctic dipolc braking model (see the pulsar catalog by Pay 



enhanded scenarios of thermal evolution, depending on 

whether the dominant process is cither the modified or 
the direct Urea. As it was shown in Schaab et al. (1997 ) 



lor et al., 1993)). Because of the highly conducting core 
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the surface may hide an even stronger magnetic field in the 
interior (Muslimov & Page, 199£). By applying the New- 
tonian scalar virial theorem Lai & Shapiro (1991) showed 
that the interior field may be as high as 10 18 G. Works on 
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the origin and evolution of the interior magnetic field show 
that at least field strengths ot_H = 10 15 to 10 17 G are pos- 



( Bisnovatyi-Kogan & Moschcnko, 1992; Bisnovatyi- 


i 1992; 


Duncan & Thompson, 1992 




Thompson & 



Kogan 

Dunca i, 1995 ) , which is in the range of the critical field 
strength for m = 2-superfluidity. 

Besides the neutrons also the protons supposedly pair 
in the interior of a neutron star. Since they probably form 
a type II-superconductor, the magnetic field is confined to 
flux tubes, whose radial dimension is given by the London 
penetration depth, A ~ 40-100 fm. The magnetic field 
inside the flux tubes is therefore of the order of H « 
$ ln(A/£)/27rA 2 - 10 14 -10 16 G, where $ ~ 2.1 x 10~ 7 
G cm 2 is the flux quantum and £ ~ 10 fm is the coherence 
length. Even if the bulk of the neutrons do not experience 
a field above the critical field H c , the field inside the flux 
tubes may be large enough. This would be sufficient for 
a considerable modification of the thermal properties of a 
neutron star. 

In this work we will study the effect of nodes of the neu- 
tron gap energy along the magnetic axis on the thermal 
evolution of neutron stars and compare the results with 
neutron star models with 3 P2-pairing in the m = 0-state 
and also with observed data. We will start with the un- 
derlying physics by discussing the energetically favoured 
pairing state in 3 P2-superfluidity and the suppression of 
the most important neutrino emission processes (Sect. ^] 
and ||). In Sect. |^ we classify our models into 16 scenar- 
ios. The outcome of the simulations for these scenarios is 
compared with observations in Sect. ^. Section ^| contains 
the conclusions and the discussion. 



1. F2-pairing 

Neutrons are believed to form 1 So-superfluid pairs in the 
density range between neutron drip and the saturation 
density of normal nuclear matter. At supernuclear densi- 
ties, i.e. in the core regime of neutron stars, the repulsive 
interaction becomes dominant and the ^o-superfluidity 
gap closes. However, the tensor and spin-orbit interaction 
become attractive at such densities and the formation of 
3 P2-neutron pairs becomes possible. The ordering parame- 
ter for 3 P2-pairing is a complex, traceless, symmetric 3x3- 
matrix A. It can be determined in vicinity of the transition 
temperature T c by minimising the Ginzbur g-Landau free- 



energy functional f2[A]. This was done by Muzikar et al 
1980|) , who found that for a magnetic field strength 



H > H C (T) = 3 x 10 17 



knT r 



1 



T 



1 MeV 

the energy gap for quasiparticle of momentum p is 

3 



G (1) 



>(P)I' 



8tt 



A 2 l=2 (T)sin 2 



(2) 



Here, fee denotes the Boltzmann-constant, 9 the angle be- 
tween the magnetic field H and quasiparticle momentum 



p, p the unit vector in direction of p, and A 2 ra=2 (T) a tem- 
perature dependent factor. For vanishing magnetic fields 
the free energy functional is minimised by the gap function 

=o(p)| 2 = r±-tf n=0 (T)(l + 3cos 2 0') , (3) 



\A r 



16tt 



where 9' is the angle between p and some quantisation 
axis. As we have already noted the gap function (|^) van- 
ishes along the magnetic axis at 6 = 0, ir, whereas the 
gap function (||) is nodeless. For magnetic fields < H < 
H C {T) the minimising gap function is a superposition of 
Eqs. ^ and |[ which is still nodeless. 

The critical field strength vanishes at the critical tem- 
perature T c and increases to 



ff c (T = 0) = 3x 10 17 ^kG 



(4) 



lMeV 

at T — 0. In the pairing model of iHoffberg et al. (1970| ) the 
maximum energy gap is, for instance, A max fa 4.5 MeV. 
This yield a critical field strength H C (T = 0) « 2 x 10 17 G. 
If the inner magnetic field is smaller than H C (T = 0), first 
a m = 2-superfluid is formed at T c and then at 

H 



Tr 2 =T C \ 1 — 



H c (T = 0) 
a nodeless m = 0-superfluid is formed. 



(5) 



3. Neutrino emissivity in superfluid matter 

In the vicinity of the Fermi surface (\p— pf\ *C pf) the one- 
article energy can be approximated by (see, e.g., Lifshitz 
Pitaevskii, 1980|) 



E(p) 



Pf) 



^\A(pW+vUp 



Pf) 



for p < p Fl 
for p > p F , 



(6) 



where up, p-p, and fi is the Fermi velocity, momentum, 
and energy, respectively. The superfluidity creates a gap 
of width 2|A(p)| in the energy spectrum around the Fermi 
momentum p-p. Excitations of particles above the Fermi 
energy are therefore suppressed for fceT < |A(p)|. We shall 
study the two possible cases of triplet pairing of neutrons 
which occur (or might occur) in neutron stars' interior (see 
Sect. ^). The gap function |A m= o(p)| in the case of node- 
less pairing (Eq. |3|) has the approximate form (Lcvenfish 
|fc Yakovlev, 1994b[ ) 



Mi» = y ( 0.7893+ 1.188^ ) L- L ;T 



x ^/l + 3co^9, (7) 



whereas gap (H) gives 



|A(p)r 



= 11- 



>, | 2.(13(1 -(UJ)()3( +0.1727 (^-j J^sm, 



(8) 
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The suppression of the neutrino emissivity e u by su- 
perfluid neutrons can be expressed by a factor 



c 

,n.sf ' 



(9) 



which was calculate d for a couple of processes in Levenfish 
& Yak( |)vlev (1994b[) and |Yakovlev fc Levenfish (1995|). Let 



us c onsider for example th e direct Urea process ( Boguta. 
198l |; l^attimcr et al., 199l| ) 



n-*-p + e +D C . (10) 

The suppression factor can analytically be calculated only 
in the limiting case T -> ( |Lcvcnfish fc Yakovlcv, 1994b| ) 



126 
457tt 5 



602971-" -2 
5484 



for case m = 0, 
for case m = 2, 



(11) 



where the dimensionless parameter v is given by 



1 A(T) 



167T fe B T 
3 A(T) 



for case m = 0, 
for case m = 2. 



(12) 



By numerical fitting Levenfish & Yakovlev (19941) obtain 



(o.2546 + ^(0.7454)2 + (0.1284i;) 2 
x exp ^2.701 - v/(2.701) 2 + v 

| + (0.09226ti) 2 
1 + (0.1821d) 2 + (0.16736d) 4 



I exp ( 1 - a/1 + (0.4129?;)' 



for m = 0, 



for m 



(13) 

As can be deduced from Eq. [ll|, the suppression factors of 
the nodeless case m — decreases with decreasing tem- 
perature much more rapidly than the suppression factor 
for case m — 2 (see Fig. [I]). This is caused by the two 
nodes of the gap function at the intersections of the Fermi 
surface with the magnetic axis. 

Besides the direct Urea process we considered a couple 
of oth er proces ses (see Table [|). Following Levenfish fc 



Yakovl ;v (1996 ) we approximate the suppression factors 



for these processes by 



(14) 



if only one neutron is involved (proton branch of mod- 
ified Urea process, electron-neutron, and neutron-proton 
scattering), and by 



n 2 n (v) = n d n (2v) 



(15) 



if two neutrons are involved (neutron branch of modified 
Urea process and neutron- neutron scattering). 

The suppression factors for proton superconductivity 
are given in Yakovlcv fc Levenfish (1995| ) . If both neutrons 
and proton are superfluid we approximate the combined 
suppression factor 72. np by the minimum of the factors lZ n 
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Fig. 1. Suppression factor TZ^ of the direct Urea process 
for the two cases of triplet pairing. 



and 1Z 

other nuclear species (cf. Levenfish 



p obtained by neglecting the superfluidity of the 

Yakovlev, 1996| ): 



mm(TZ n ,TZ p ). 



(16) 



We account also for the suppression of the heat capacity 
( Levenfish fc Yakovlev, 1994a) an d of the heat conduc- 
tivity (|Gnedin fc Yakovlev, 1995|) . Thought the factors 
are similar as in the case of neutrino processes, the effect 
on the thermal history is only small compared to the ef- 
fect of the suppression of the neutrino emissivity. In the 
case of the heat capacity at least electrons and possibly 
also non superfluid protons and hyperons still contribute 
to the total heat capacity, whereas the heat conductiv- 
ity is maintained by electron-electron, electron-muon, and 
muon-muon scattering processes. 



4. Cooling scenarios 

The general relativistic equations of stellar structure 
and thermal evolution (cf. Thornc, 1977) were numeri- 
cally solved via an implicit finite difference scheme by a 



Newton- Raphson algorithm (see Schaab et al., 1996 for 



more details). The outcome of the simulations depend on 
a couple of ingredients which are partly rather uncertain. 
If one took all these uncertainties into account, one would 
get a huge number of models, which are difficult to sur- 
vey. However all these models can be divided into four 
disti nct classes dependin g on the dominant neutrino pro- 
cess ( [Schaab et al., 19*97 ). These are: 



Class A - Standard cooling: This class consists of slowly 
cooling models which do not allow for any enhanced 
cooling process (like the direct Urea process). 

Class B - Intermediate cooling: The emission rate of the 
modified Urea process is increased by medium effects 



( Voskresenskii & Senatorov, 1986; 3chaab et al., 1997) 
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Table 1. Neutrino processes considered in the cooling simulations 



Process 



Emissivity from Ref. 



pair-, photon-, plasma-processes 
bremsstrahlung in the crust 
bremsstrahlung in the core: 

n-n scattering 

n-p scattering 

p-p scattering 

e~-n scattering 

e~-p scattering 

e~-F scattering 
superfluid pair breaking and formation 
modified nucleon Urea: 



[toh et al (1989D 



Yakovlev fe Kaminker (1996); Haensel et al. (1996) 



Friman k Maxwell (197S ) 



Friman k Maxwell (197S ) 



Friman k Maxwell (1979T ) 



Kaminker et al. 



Kaminker et al. 



(1997 ) 



Kaminker et al. 



(1997 ) 



(1997D 



Voskresenskii k Senatorov (1987); schaab et al. (1997); f 



neutron branch 


Friman k Maxwell (1979) 


or alternatively 


Voskresenskii k Senatorov (1986); schaab et al. (1997) 


proton branch 


Vakovlev k Levenfish (19951) 


direct nucleon Urea 


Lattimer et al. (1991) 


direct hyperon Urea 


Prakash et al. (1992|) 



t: In the case of P2-pairing the emissivity was recalculated by averaging over the angle 9 between quasiparticle 
momentum p and quantisation axis (s. Eqs. ^| |^) 



This leads to a considerably faster cooling during the 
neutrino cooling era. 
Class C - Enhanced nucleon cooling: The third class al- 
lows for the direct nucleon Urea process (maybe in 
combination with pion or kaon condensation) yielding 
an even faster cooling than in the intermediate cooling 
scenario. 

Class D - Enhanced hyperon cooling: The last class al- 
lows additionally or exclusively for the direct hyperon 
Urea process, which is not suppressed by nucleon su- 
perfluidity in contrast to the nucleon processes^} 

The possibility of the direct nucleon and hyperon Urea 
process depends on the composition in the interior of the 
star and hence on the underlying equation of state and 
on the central density of the star. The equations of state 
can therefore also be divided into these four classes, if the 
mass of the neutron star model is fixed to some value (we 
selected the canonical value M = 1.4M Q ). Such a division 
is done in Table |[ We used a broad collection of equa- 
tions of state which comprises rclativistic field-theoretical 



the effective pion gap - which are unfortunately rather 
uncertain. 

As already mentioned in Sect. ^ the second important 
uncertainty concerns the superfluidity of neutrons in the 
interior of the star. Again we classify the possible models 
into four classes: 

Class 1: Neutrons pair in the nodeless 3 P2-state (m = 0). 

Class 2: The inner, uniform magnetic field strength H ex- 
ceeds its critical value H C (T = 0) « 2 x 10 17 G (see Eq. 
[j]) even for zero temperature. The neutrons therefore 
form a m = 2-superfluid. 

Class 3: The inner magnetic field is set to H = 10 17 G < 
H c . The critical field strength is therefore exceeded for 
temperatures T c .2 <T<T C only. 

Class 4: Either the critical temperature of neutron pairing 



is very small (e.g. Elgar0y et al., 1996a), or the super- 
fluid phase does not extend over the whole interior of 



the neutron star (e.g. Takatsuka, 1972) 



)dcL 



one non-relativistic Thomas-Fermi model, as well 



For the classes 1-3 we used the gap energy by Hoffberg 
et al. (197C )^J. For class 4 we used the models by Elgar0y 



as nonl-relativistic Schrodinger-ha,sed ones. The first two 
classes, the standard and the intermediate cooling classes, 
base on the same set of equations of state, but on differ- 
ent neutrino emission rates of the modified Urea process. 
The division in two different classes seems sensible, since 
the rate of this process depends quite sensitively on the 
choice of the microscopic parameters - as, for instance, 



et al. (1996a) and Takatsuka (1972). The angle-averaged 
energy gaps of 3 P2-pairing in the m = 0- and m = 2- 
states are almost equal (Elgar0y et al., 1996t). We can 
therfore approximate the energy gap for the m = 2- 
state, which is not given by all authors, by A m= 2(T) « 



iA m=0 ( Amundsen fe Cstgaard, 1985 ) 



The ability of hyperons to form superfluid or 



2 Hoffberg et al. (197C) estimated the gap energy of 3 P2- 



ducting pairs has only recently been investigated (Balberg k 
Barnea, 1998). We shall assume here that hyperons cannot pair 



supercon- 



and a ddress this problem in a forthcoming work (Schaab et al. 
199Sl). 



pairing by using the phase shift data of neutron-neutron scat- 
tering and assuming the effective nucleon mass to be equal 
to its bare mass. The obtained gap energy serves therefore as 
an upper limit and is used here to demonstrate the effect of 
rn — 2-pairing. 
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Table 2. Classification of equations of state into the four cooling classes denned in the text. The canonical neutron 
star mass M = 1.4M© is assumed. The classification into the standard and intermediate cooling class depends on 
whether or not medium modification of the modified Urea process are taken into account. 



Class 


JiUo 


Baryonic 
composition 




Interaction 


Many-body 
approach 


Ket. 


A/B 


UV14+UVII 


p,n 




2 nuclei potential Urbana V14 and 
3 nuclei potential Urbana VII 


NRV 


1 


A/B 


UV14 + TNI 


p.n 




2 nuclei potential Urbana V14 and 
density dependent terms TNI 


NRV 


1 


A/B 


TF 


p,n 




2 nuclei potential TF 


NRTF 


2 


C 


HV pn 


p,n 




exchange of <J,u,p mesons 


RH 


3 


C 


RHFp n 


p,n 




exchange of a, lo, p, it mesons 


RHF 


4 


c 


/~iK240 
^pn 


p,n 




exchange of o,uj,p mesons 


RH 


5 


D 


HV 


P,n,A I E ± . u J H". 




exchange of a, u),p mesons 


RH 


3 


D 


RHF1 


p^AE^ ,- -- 


•A 


exchange of a,u,p,ir mesons 


RHF 


4 


D 


RHF8 


p.n.A.E^ ^" 


A 


exchange of a,u,p,TT mesons 


RHF 


4 


D 


/-1K240 
^^478 


p,n ) A ) E ± .° J H ' 




exchange of a,u),p mesons 


RH 


5 


D 


n K240 
«M78u 


P^A.E^ ,- ' 




exchange of o,uj,p mesons 


RH 


5 



Abbreviations: NRV: non-relativistic variational method, NRTF: non-relativistic Thomas-F ermi model, RH: rela - 
ti yistic Hartree approx ima tion, RHF: relativistic H art ree-Fock approxima tion. References: 1: Wiringa et al. (1988| ), 
2: |Strobel et al. (1997|), 3: |Weber fc Weigel (1989|), 4: |Huber et al. (1997|), 5: Glendenning, priv. comm. 



The cooling simulations are performed by running the 
numerical code for all combinations of parameters which 
are in agreement with the respective cooling scenario (a 
detailed list of all parameters can be found on the Web: 



http: / /www. physik.uni-mucnchcn.de/scktion/suessmann 
/astro/cool/schaab.0797/input.html). This procedure 
leads to cooling bands rather than individual cooling 
tracks, by means of which we demonstrate the uncertainty 
inherent in the input parameters. By combining the four 
classes of equations of state A-D with the four classes of 
superfluidity 1-4 we get a collection of sixteen scenarios 
whose behaviour will be studied in the next section. 

We shall study the case where the magnetic field ex- 
ceeds H C (T — 0) only inside the flux tubes of the proton 
type II-superconductor separately. In this case (labeled as 
class 27), the suppression factor 1Z can be expressed by 
the combination 

n = jiiin + (1 - 7 )7e out (17) 

of the suppression factors !Zi n , T^cmt inside and outside the 
tubes. 7 is the volume fraction occupied by the proton 
flux tubes. Inside the tubes we assume that the magnetic 
field exceeds the critical field strength H C (T — 0). Because 
protons are unpaired inside, the suppression factor 7ti n 
is equal to the suppression factor for neutron m = 2- 
pairing. Since the magnetic field is constrained to the flux 
tubes, neutron pair in a nodeless m = 0-state outside the 
tubes. lZ ou t accounts for both superfluids, the nodeless 
neutron superfluid and the singlet proton supcrfluid by the 
approximate expression of Eq. (|l^). The volume fraction 
7 is determined by the relation 

Hcx (18) 



Htu 



be 



of the bulk field strength 7J 0X and the field strength inside 
the proton flux tubes -fftube, which we set equal to the 
critical field strength H C (T = 0) w 2 x 10 17 G. 

5. Results and comparison with observations 

Among the soft X-ray observations of the 27 sources which 
were identified as pulsars, the ROSAT and ASCA observa- 
tions of PSRs 0002+62, 0833-45 (Velah 0656+14, 0630+18 
(Geminga) and 1055-52 (see Table y) achieved a suffi- 
ciently high photon flux such that the effective surface 
temperatures of these pulsars co uld be extracte d by two- 
or three-component spectral fits ( Ogelman, 1995 ). The ob- 
tained effective surface temperatures, shown in Figs. || to 
[| depend crucially on whether a hydrogen atmosphere is 
used or not. The kind of atmosphere of individual pulsars 
could be determined by considering multi-wavelength ob- 
servations ( Pavlov et al., 1996 ). Since the measured spec- 
tra are however restricted to the X-ray energy band up to 
now, one is not in the position to determine the kind of at- 
mosphere. For that reason we investigated both the black- 
body model and the hydrogen-atmosphere model, distin- 
guished in Figs. ^|to|^ by error bars with a solid and hollow 
circle, respectively. All error bars represent the 2a error 
range due to the small photon fluxes. 

Except for PSRs 0833-45 (Vela) and 0002+62, all ages 
are estimated by their spin-down age r = P/2P. This 
relation implies however that both the moment of inertia 
and the magnetic surface field are constant with time, and 
that the braking index n is equal to its canonical value 
3 (angular momentum loss due to pure magnetic dipole 
radiation) . The true ages may therefore be quite different 
from the spin-down ages. The age of Vela was recently 
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Table 3. Surface temperatures as measured by an observer at infinity, T s °°, and spin-down ages, r, of observed pulsars. 



Pulsar logr [yrs] 



Model atmosphere log T s °° [K] 



Reference 



0002+62 


~ 4t 


blackbody 6.20^;^ 


Hailey & Craig (1995), Fig. 2 with lower limit on jV H 


0833-45 


4.4 ±0. If 


blackbody 6.24 ± 0.03 


Ogelman (1995), Table III 



(Vela) 



magnetic H-atmosphere 5.85 ± 0.09 



0656+14 5.05 blackbody 5.89 

magnetic H-atmosphere 5.72 



+o.o« 

-0.17 
+0.06 

-0.05 



r 7 r+U.Ub 
°- l,J -0.08 

5 42+ ' 12 

°- 4Z -0.04 



Greiveldinger et al. (1996;), T able 2 
lAnderson et al. (1993;) 



|Halpern fc Wang (1997] ), Table"! 
Meyer et al. ( 1994j), Fig. 2a with B 12 = 1.18 



0630+18 5.53 blackbody 

(Geminga) magnetic H-atmosphere 



1055-52 



5.73 



blackbody 



5.90loT 2 t 



Greiveldinger et al. (1996), Table 2 



f estimated true age instead of spin-down age (see text). 

J The 2cr range is not given in the respective reference. We therefore estimated it by multiplying the given la 
error with a factor 2. 



determined by Lync et al. (1996 ), and the approximate 
age of PSR 0002+62 is given by an estimate o f the age 
of the relate d supernova remnant G 117.7+06 ( Hailey fc 
Craig, |1995|). 



All calculations were performed for a gravitational 
mass of M — 1.4M , about which the observed pulsar 
masses tend to scatter. In Figs. || to rej, we plot the surface 
temperature of the neutron star models, as observed at 
infinity, against the stars age. Figure || shows the cooling 
bands for superfluidity case 1 and the equation of state 
classes A-D (see Sect. || for the definition of the classes). 
Superfluidity suppresses the nucleon processes so strongly, 
that the bands of the three classes standard, intermediate 
and enhanced nucleon cooling nearly coincide. Since the 
hyperon processes are not suppressed within the models 
used here, the band of scenario Dl is located much be- 
low the other three bands. Note that the relatively low 
surface temperature, in comparison with other works (e.g. 
Umcda et al., 1994| ; [Page, 1995| ; |Schaab et al., 1996| ), of 
the standard scenario Al is caused by the inclusion of the 



superfluid pair breaking and formation process (Voskre 



senskii & Senatorov, 1987; Schaab et al., 1997), which is 



not included in older works. The surface temperature ob- 
tained by fitting with a magnetic hydrogen atmosphere 
(hollow circles) lie right between the scenario Dl and the 
scenarios Al-Cl. In this model of superfluidity the gap 
between the bands can only be filled by redu cing the su- 
perfluid critical tem perature of the neutrons ( Page, 1995 ; 



Schaab et al., 1996[ ) 



By increasing the inner magnetic field strength above 
its critical value we reach the second case of superfluidity 
(Fig. ||). Band C2, which corresponds to enhanced nucleon 
cooling, moves to lower temperatures. The observation of 
PSR 0630+18 (Geminga) is now consistent with this sce- 
nario as long as this pulsars has a magnetised hydrogen 
atmosphere. 

If the magnetic field is in the order of the critical field 
strength but still below it, am = 2-superfluid is first 
formed by the neutrons, and later (as the star cools below 




3.0 4.0 

log(T/yr) 

Fig. 2. Thermal evolution of a M — 1.4M Q neutron star 
with superfluidity class 1. The bands correspond to the 
four classes of the equation of states as described in Sect. 



?c, 2, see Eq. [|) this superfluid turns into a nodeless state. 
This possibility is referred to by class 3 (Fig. ||). Only in 
scenarios C3 and D3 the temperature reaches this second 
critical value T C 2 during the neutrino cooling era, that 
is during the time where cooling is dominated by energy 
loss due to neutrino emission. For that reason the band 
of the C3-scenario is pushed up again to higher temper- 
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log(x/yr) 

Fig. 3. Same as Fig. || but for superfluidity class 2. 

ature which are consistent with the surface temperature 
obtained by fitting the observations with a magnetised hy- 
drogen atmosphere (apart from PSR 0630+18). Scenarios 
B3 and A3 show again the same cooling behaviour. 

In case 4, we assume that the neutron superfluid does 
not extend throughout the whole neutron star's core or 
the critical temperature is very small T c < 10 s K (Fig. 
^). Then the nucleon processes are not suppressed every- 
where, and particularly the models of the C4-scenario of 
enhanced nucleon cooling cool down very fast. Certainly 
this scenario, as well as scenarios D1-D4, is not consistent 
with the observations, no matter what kind of atmosphere 
is used. Case 4 would therefore clearly favour the standard 
or the intermediate cooling scenarios. The higher emissiv- 
ity of the modified Urea process in scenario B4 lead now to 
considerably smaller surface temperatures in comparison 
with scenario A4. 

Finally we consider the case 27, where the critical field 
strength H C (T = 0) is exceeded only inside the flux tubes. 
The suppression factor TZ in Eq. (|lj]) depends on the vol- 
ume fraction 7 occupied by the flux tubes. We study the 
cooling behaviour for two values of 7: 7 = 10~ 4 , which cor- 
responds to a bulk magnetic field strength H — 2 x 10 13 G 
(see Eq. |l8|), and 7 = 10~ 2 , which corresponds to H = 
2 x 10 15 G. To clarify the impact of this case on the cool- 
ing history of neutron stars, we plot cooling tracks for the 




log(i/yr) 

Fig. 4. Same as Fig. || but for superfluidity class 3. 

HV pn equation of state instead of bands (Fig. ||). Since the 
neutrino emission is exponentially suppressed only outside 
the proton flux tubes the effect of suppression is not as 
high as in the case of model CI (which corresponds to 
7 = 0), but still larger than in the model C2 (which cor- 
responds to 7 = 1). The obtained surface temperatures 
are consistent with the surface temperature obtained by 
fitting the observations with a magnetic hydrogen atmo- 
sphere (again apart from PSR 0630+18). Although the 
surface temperature obtained in both scenarios C3 and 
C27 with 7 = 10~ 2 are consistent with the observation of 
PSR 0656+14 within the magnetic hydrogen atmosphere 
model, the bulk magnetic field strength in both models 
differ by two orders of magnitude. 

6. Conclusions and discussion 

In this work we have incorporated the possibility of nodes 
of the neutron gap energy along the magnetic axis in the 
core of a neutron star in simulations of thermal evolution. 
The realization of such superfluids depends on a rather 
strong magnetic field, that might however occur in the 
interior of a neutron star or in the interior of flux tubes of 
the proton type II superconductor. In any case the state 
of m = 2-superfiuid turned out to have several interesting 
consequences. 
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3.0 4.0 

log(t/yr) 

Fig. 5. Same as Fig. but for superfluidity class 4. 



Since the suppression factor 1Z of neutrino emissivity 
decreases only polynomially as function of T/T c instead 
of exponentially, the effect of superfluidity is not as large 
as in the case of nodeless (m = 0) neutron superfluidity, 
which was studied in literature so far. Particularly the 
differences in the undiminished neutrino emissivities of the 
cooling classes A-D remain visible in the resulting surface 
temperatures. 

The second important consequence is that in scenario 
C3 the second critical temperature T c ^ depends on the in- 
ternal magnetic field (see. Eq. [""]). In scenario C27 the sup- 
pression factor depends on the factor 7 (see Eq. |l7| ), which 
itself depends on the bulk magnetic field strength (see Eq. 
|l8| ) . By varying the magnetic field the band corresponding 
to scenarios C2j and C3 can be moved to lower or higher 
surface temperatures in the age-range 10 2 yr < r < 10 5 yr. 
This can also be achieved by varying the maximum en- 
ergy gap, as it was done by Page (1995 ) and Schaab et al 
(1996[)] Provided, that the energy gap reaches its maxi- 
mum value for densities smaller than the central density 
of a typical neutron star, the maximum energy gap should 
be the same for all neutron stars. This is different with re- 
spect to the inner magnetic field strength since at least 
the surface magnetic field strengths of individual pulsars 
differ by five orders of magnitudes. This consequence may 
become important if two neutron stars are found which 



3.0 4.0 

log(x/yr) 

Fig. 6. Thermal evolution of a M = IAMq neutron star 
obtained for the HV pn equation of state. The curves cor- 
respond to different superfluidity classes. 



have about the same age but different surface tempera- 
tures, thought this might also be explained by different 
star masses in connection with the intermediate cooling 
scenario (see Fig. 8 in fjchaab et al., 19*97 ) and by more 
or less accreted matter on the surface of the star (see Fig. 
9 in Potckhin et al., 1997). The two scenarios C27 and 



C3 differ in the internal magnetic field strength which is 
needed to explain, e.g., the surface temperature of PSR 
0656+14 obtained within the hydrogen atmosphere model. 
Whereas for scenario C3 the rather strong magnetic field 
H ~ 10 17 G is needed, a bulk magnetic field strength 
H r~ 10 15 G is sufficient in scenario C27. 

The possibility of the direct nucleon Urea process 
becomes more and more probable. To the best of our 
knowledge all relativistic Hartree-, Hartree-Fock-, and 
Bruckncr-Hartree-Fock-equations of state have sufficiently 
large proton fractions to allow for the direct Urea pro- 
cess. In the past there were some doubt about the relia- 
bility of such high proton fractions, because the nonrel- 
ativistic variation al equations of state, like UV14+UVII 
and UV14+TNI ( |Wiringa et al., 1988j ), predict proton 
fractions which are below the threshold value for the di- 
rect Urea process, since the symmetry energy does not 
monotonously increase with density. This seems to depend 
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however on the used potentials, since more modern po- 
tentials, like AVi8, Nijm-1,11, Rcid93 and CD-Bonn, yield 



higher proton fractions ( Engvik et al., 1997 ). 

It seems therefore reasonable to restrict our following 
discussion to classes C and D. Since the impact of hyperon 
pairing has not been investigated yet, let us assume that 
all observed neutron stars belong to class C. Under this as- 
sumption the situation would become rather enlightened 
if the chemical composition of the neutron stars envelope 
could be determined by multi-wavelength observations. In 
the case that the blackbody fits will turn out to be ade- 
quate only scenario CI is in agreement with the observa- 
tions. If the magnetic hydrogen fits are more reasonable 
only scenarios C2, C27, and C3 are consistent with the 
data. However also scenarios CI and C4 can be brought 
into agreement i f one reduces the ga p energy in the case of 
CI ( Page, 1995 ; Bchaab et al., 1996 ) or takes eff ects of an 



accret ed en v elope into a ccount in the case of C4 ( [Potckhin 
et al., jHf ; |Page, 1996|) . As it was already claimed, these 



different scenarios could be further distinguished, if the 
surface temperature of two or more neutron stars with 
similar ages could be determined. 
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